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Abstract 

We re-examine Peres' statement "opposite momenta lead to opposite 
directions" . It will be shown that Peres' statement is only valid in the large 
distance or large time limit. In the short distance or short time limit an 
additional deviation from perfect alignment occurs due to the uncertainty 
of the location of the source. This error contribution plays a major role in 
Popper's orginal experimental proposal. Peres' statement applies rather 
to the phenomenon of optical imaging, which was regarded by him as a 
verification of his statement. This is because this experiment can in a 
certain sense be seen as occurring in the large distance limit. We will also 
reconsider both experiments from the viewpoint of Bohmian mechanics. 
In Bohmian mechanics particles with exactly opposite momenta will move 
in opposite directions. In addition it will prove particularly usefull to use 
Bohmian mechanics because the Bohmian trajectories coincide with the 
conceptual trajectories drawn by Pittman et al. In this way Bohmian 
mechanics provides a theoretical basis for these conceptual trajectories. 

1 Introduction 

Already in 1934 Popper [1-3] proposed an experiment which aimed to test the 
general validity of quantum mechanics. Popper assumes a source S from which 
pairs of particles are emitted in opposite directions. Two observers, say Alice 
and Bob, are located at opposite sides of the source, both equiped with an array 
of detectors. If Alice puts a screen with a slit in her way of the particles, she will 
observe a diffraction pattern behind the screen. According to Popper, quantum 
mechanics will also predict a diffraction pattern on the other side of the source, 
where Bob is located, when coincidence counts are considered. This is because 
every measurement by Alice is in fact a virtual position measurement of the 
correlated particle on Bob's side, leading to an increased momentum uncertainty 
for Bob's particle as well. This is the same diffraction pattern that would be 
observed when a physical slit was placed on Bob's side. Popper, who declared 
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himself a metaphysical realist, found this idea of "virtual scattering" absurd 
and predicted no increased momentum uncertainty for Bob's measurement due 
to Alice's position measurement. He therefore saw his proposed experiment as 
a possible test against quantum mechanics and in favour of his realist vision in 
which particles have at each time well defined positions and momenta; particles 
for which the Heisenberg uncertainty, for example, is only a lower, statistical 
limit of scatter. 

Unfortunately, to describe the setup of the experiment, Popper occasionally 
invoked classical language, which veiled some severe problems which could ob- 
struct a practical realization of his experiment (for an extensive discussion see 
Peres [4]). Popper writes: "We have a source S (positronium, say) from which 
pairs of particles that have interacted are emitted in opposite directions. We 
consider pairs of particles that move in opposite directions . . . " . It is the valid- 
ity of this statement, which appears to be a very delicate issue, that we shall 
deal with in the present paper. 

Of course, when we consider a decaying system at rest in classical mechanics, 
the fragments will have opposite momenta 

Pl+P2=0 (1) 

and if we take the place of decay of the system as the centre of our coordinate 
system the positions of the two fragments will satisfy miri + m 2 r 2 = 0, so the 
fragments will be found in opposite, isotropically distributed directions. If the 
fragments have equal masses then they will be found at opposite places, relative 
to the centre of the coordinate system. 

But these properties don't hold in quantum mechanics. Suppose that a 
system has a two-particle wavefunction ip which has a sharp distribution at 
P1+P2 = 0, i.e. both \(ftij+P2j)\ and A(pij +P2j) are small for every component 
j of the momentum vectors pi +p2- Then according to the uncertainty relations 

A(pij +p 2 j)A(m 1 fij +m 2 r 2 j) > +m 2 ), (2) 

the distribution of miry + m 2 r 2 j will be broad for every component j. In the 
case of equal masses, the inequalities in Eq. (2) imply that opposite momenta 
are incompatible with opposite positions. In particular, at the moment of decay, 
the inequalities imply that opposite momenta of the particles are incompatible 
with the latter being both located at the origin of our coordinate system. It 
was even shown by Collett and Loudon [5] that this initial uncertainty on the 
location of the source implies that Popper's experiment is inconclusive. 

Although the original proposal of Popper's experiment can hence not be 
performed practically, due to the fact that opposite momenta are incompatible 
with opposite positions, the intention of Popper's proposal can be maintained 
if we have a two-particle system which displays some form of entanglement in 
the position coordinates. 1 One is then in principle able to test experimentally 

1 The position entanglement should not be exact, otherwise, as shown by Short [6], both of 
the observers would observe an infinite momentum spread. 
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wether one of the particles will experience an increased momentum spread due 
to a position measurement (within a slit width) of the correlated particle, i.e. 
we would be able to test a possible "virtual scattering". Such an experiment 
has actually been carried out by Kim and Shih [7], who used the phenomenon 
of optical imaging, which was reported by Pittman et al. [8] . In the process of 
optical imaging, position entangled photons are created by placing a lens in the 
way of one of the momentum entangled photons created by parametric down 
conversion. 

Recently Peres gave an analysis to which extent opposite momenta lead 
to opposite directions [9]. He argues that the inequalities in Eq. (2) do not 
exclude a priori the possibility that opposite momenta of particles lead to op- 
posite directions (instead of positions) where the particles will be found when 
a measurement is performed. On the contrary, the operator equivalent of (1) 
would even lead to an observable alignment of the detection points of the two 
particles. However in section 2 we show that Peres' analysis can only be ap- 
plied in what we will call the "large time" or "large distance" regime. This is 
the regime where the two particles have travelled a "large" distance from the 
source. In this limit the Heisenberg uncertainty will be of minor importance for 
the angular correlation. This result agrees with the scattering into cones theo- 
rem, which applies in the limit of infinite time. In the "short" distance regime 
the uncertainty on the location of the source, due to the Heisenberg uncertainty, 
leads to an additional error contribution to the angular alignment. We will also 
make clear what is exactly meant by the "large" or "short" distance regime. 

An example which shows the importance of the additional error, even in the 
long time regime, is the orginial experiment of Popper (section 3). The situation 
is different for the phenomenon of optical imaging, which was considered by 
Peres as a justification for his analysis. Because the experiment can in a certain 
sense be seen as occuring in the long distance limit, Peres' analysis can be 
applied to this experiment. 

In section 4, we will provide a simplified description of a decaying system us- 
ing the Bohmian picture [10, 11]. It will prove usefull to use Bohmian mechanics 
because Bohmian mechanics inherits some of the language of classical mechanics. 
This is because Bohmian mechanics describes the motion of particles, allowing 
one to consider possible trajectories of the particles. This is contrary to "stan- 
dard quantum mechanics" where we can only speak of probabilities where a 
particle will be detected when a position measurement is performed. Before the 
position measurement is performed, particles can't, according to standard quan- 
tum mechanics, be seen as localized pointlike objects. As such, there is no basis 
in the standard quantum formalism to use trajectories of particles between two 
measurements. The classical language present in Bohmian mechanics can then 
be compared with the classical language used by Popper in his experimental 
proposal, in order to examine the possible interplay between opposite momenta 
and opposite positions. 

We will also give an alternative explanation of the phenomenon of optical 
imaging (or at least of its massive particle equivalent) with the use of Bohmian 
mechanics (section 5). Although the experiment could be equally well described 
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by quantum optics, there is an additional advantage by using Bohmian mechan- 
ics. The paper of Pittman et al. contains drawings of conceptual trajectories 
of the photons. These trajectories are not the real paths followed by the pho- 
tons, as explained above, but merely serve as a tool to visualize the experiment. 
In section 4, we show that the Bohmian trajectories (calculated in the massive 
particle equivalent of the experiment of Pittman et al), coincide with these con- 
ceptual trajectories in Pittman's paper. As such, Bohmian mechanics provides 
a theoretical basis for these trajectories. 

Wc want to stress that by using Bohmian mechanics we don't intend to 
contest the validity of standard quantum mechanics. Besides, as is well known 
Bohmian mechanics and standard quantum mechanics are completely equivalent 
at the empirical level. Now, wether one believes or not in the existence of 
particles as localized objects between two measurements, Bohmian mechanics 
provides a possible way of dealing with particle trajectories, consistent with 
quantum mechanics at the level of detections. Even if one is dissatisfied with 
the picture of moving particles, one can interpret the Bohmian trajectories as 
the flowlines of the probability, because the speeds of the Bohmian particles are 
defined as proportional to the quantum mechanical probability current. 

Finally we want to remark that, unfortunately Kim and Shih failed in their 
original intention to perform Popper's gedankenexperiment. It was shown by 
Short [6], that due to imperfect momentum entanglement of the parametric 
down converted photons (caused by a restricted diameter of the source), the 
image was not perfect, which blurred the predicted results. 

2 Opposite momenta and opposite directions 

To discuss the possible angular alignment of momentum entangled particles, 
Peres considers a nonrelativistic wavefunction describing massive particles. Ac- 
cording to Peres, the reason for angular correlation of momentum entangled 
photons (as in the experiment of Pittman et al.) is the same as in the consid- 
ered massive case. 

The momentum correlated particles can be assumed to result from a decaying 
system at rest. The decaying system can then be described by the wavefunction 

V(ri, r 2 , t) = J F( Pl ,p 2 y^+^- Et ^ h d Pl dp 2 (3) 

where the momentum distribution F is peaked around pi + P2 =0 and around 
the rest energy Eq of the decaying system. According to Peres the opposite 
momenta of the particles lead to opposite directions where the particles will 
be found when a measurement is performed. His argument goes as follows. 
The main contribution in the integral in Eq. (3) comes from values pi and p 2 
for which p! + p 2 ~ 0. Because of the rapid oscillations of the phase in the 
integrand in Eq. (3), the integral will be appreciably different from zero only if 
the phase is stationary with respect to the six integration variables pi and p 2 
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in the vicinity of Pi + P2 = 0, i.e. 

dS dE . „ „ 

dpi dp t 

where S is the phase of F(pi, p 2 ) measured in units K, i.e. 

F(p 1 ,p 2 ) = \F(p u p 2 )\e iS ^» h (5) 

and the equations (4) have to be evaluated for Pi + P2 = 0. The equations (4) 
then determine the conditions on in order to have a non-zero ip (and |'0| 2 )- 

Peres then introduces spherical coordinates to describe Pi and , and varies 
the phase with respect to the six spherical variables Pi,(fri,6i (the spherical 
coordinates of pj). Peres further assumes that the phase of F obeys 

ds/dp l = (6) 

in the vicinity of pi + P2 = 0. This would restrict the place of decay near the 
origin of the coordinate system, because of Eq. (4). By varying with respect to 
the momentum angles, Peres obtains that the phase is stationary if p, and rj 
have the same direction. Because F is peaked at pi + P2 = 0, this results in 

9[+6' 2 = TT 

\4>'i-<&\ = tt (7) 

where r,,^,^ are the spherical coordinates of rj. These equations show that 
the two particles can only be detected at opposite directions relative to the 
centre of our coordinate system. This is because the integral in Eq. (3) (and 
hence IV'j 2 ) would only be appreciably different from zero if obey Eq. (7). 

Peres mentions two error contributions to the angular alignment. The first 
is a transversal deviation of the order yjht/m due to the spreading of the wave- 
function, which was recognized as the standard quantum limit [12]. The second 
is an angular spread of the order A(pij+p2j)/pi. One of the aims of the present 
paper is now to show that there is another error contribition which arises from 
the uncertainty on the source and which is particularly important in the "short" 
distance or "short" time regime. 

First we note that there is, apart from the conditions on 9[ and ^>\, also a 
condition on the variables ri, which is not mentioned by Peres. This condition 
is obtained by varying the phase of the integrand with respect to pi, having in 
mind the previous result that p, and r 4 have the same direction. If we define 
Vi = dE/dpi, then the additional condition reads 

n = v t t (8) 

where the Vi have to be evaluated for pi + P2 = 0. Because ip obeys the non- 

2 2 

relativistic Schrodinger equation, E equals ^7 + 2m7' w ^ n m i the masses of 
the particles. In this way Eq. (8) becomes 

n = — t, r 2 = — t. (9) 
mi m 2 
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Using Pi + P2 = one obtains 

nmi = r 2 m 2 . (10) 

Note that we haven't yet used the fact that F is peaked around a certain energy 
Eg, as is required in the case of a decaying system at rest. As Peres remarks in 
his paper, a restriction of the energy to Eq further restricts the momenta of the 
particles to satisfy pi = p\ = 2Eomim 2 /(mi + m 2 ). 

Combining (7) and (10) one obtains that the joint detection probability has 
a maximum for the classically expected relation m\V\ +m 2 r 2 = 0. Hence Peres' 
statement "opposite momenta lead to opposite directions" may be replaced by 
a stronger statement, i.e. the opposite momenta lead to a maximum detection 
probability for m\T\ + m 2 r 2 = 0. However we still have to consider the possible 
sources of deviation from this classical relation. Classically one can, in theory, 
make both quantities A(pij + p 2 j) and A(mifij + m 2 f 2 j) as small as wanted. 
Quantum mechanically one can at best prepare the system, such that initially 
the equality in 

A(pij + p 2 j)A(m 1 rij + m 2 f 2j ) > ^("M + m 2 ) (11) 

is reached. Remark that this equation implies that in case of opposite momenta, 
the particles can't depart from a confined, fixed source. 

Because the operator p\j + p 2 j commutes with the free Hamiltonian, the 
variance of the momentum operator pij + p 2 j is stationary. The variance of 
fij + f 2 j however, will in general increase with time due to the spreading of 
the wavefunction. This can be seen if we write down the expression for the free 
evolution of the operator m\Y i + m 2 f 2 in the Heisenberg picture 

mifi(t) + m 2 f 2 (t) = mifi(0) + pi(0)t + m 2 f 2 (0) + p 2 (0)i (12) 

The variance of this operator for an arbitrary component j is 

A(mifij(i) + m 2 f 2 j(t)) = A(mihj(0) + m 2 r 2 j (0)Y + A(p l3 (0)+p 23 (0))t 2 

+ (\ (mify(0) +m 2 r 2 3 

-2(mif y (0) +m 2 f 2j (0))(p lj (Q) +p 2j (0))t (13) 

where the brackets {, } denote the anticommutator. If we assume a distribution 
F which is real and symmetric, i.e. F(pi, p 2 ) = F(— pi, — p 2 ) then the last two 
terms in Eq. (13) are each zero. So the variance of m\fij{t) + m 2 r 2 j{t) increases 
with time 

A(mif y (t) +m 2 r 2j {t)) 2 = A(m 1 f lj (0) + m 2 f 2j {0)) 2 + A(p lj {0)+p 2j (0)) 2 t 2 . 

(14) 

This leads to an increasing deviation from the relation rniri + m 2 r 2 = 0. 

There is thus always an interplay between opposite momenta and opposite 
directions which is expressed in Eq. (11) and Eq. (14). We argue that one 
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should study Eq. (14), where the variances at t = in the right hand side of the 
expression are limited by the Heisenberg uncertainty in Eq. (11), to determine 
to which extent we can speak of possible angular alignment. 

Let us now see how Peres' error contributions come about and in which 
regime they are important. Assume for convenience that mi = m 2 = m. The 
transversal deviation L(t) can be taken of the order A(fij(t) + T2j(t)). There 
are two contributions to this transversal deviation. The first is A(fij(0) + 
^2j(0)) = L(0) and is important for small times. The second contribution 
is A(pij + P2j)t/m which becomes important for larger times. The angular 
deviation 9 may be derived from tan((9) = L(t)/R(t), where R(t) = pt/m is the 
distance that both particles have travelled. For small times one has tan(#) ~ 
A(rij(0) + f2j(0))m/pt and for large times one has tan(0) ~ A(pij + P2j)/p- 
Hence for large times we obtain the error mentioned by Peres. From the relations 
(11) and (14) one can also easily derive the standard quantum limit 

A(fy(t)+f 2j -(t)) 2 > 2A(fi j (0)+f 2j -(0))A(p y (0)+fty(0))t/m 

> 2ht/m. (15) 

However, this uncertainty is misleading because for small times it neglects the 
contribution arising from the uncertainty on the source A{rn\fij (0) +m 2 f 2 j (0)) . 
Especially in the considered case this contribution will be large because A (pij(0)+ 
p 2 j(0)) is small. 

In summary we see that Peres gave error contributions which only apply in 
the large time regime. These error contributions are in perfect agreement with 
the "scattering into cones" theorem which states that for every cone C in R m 
with apex in the origin 

lim / d m x\ip{x,t)\ 2 = ( d m p\(t)(p)\ 2 (16) 

with (j>{p) the momentum wave function [13]. This means that the probability 
that in the infinite future the particles will be found in the cone C is equal 
to the probability that their momenta lie in the same cone. In the small 
time regime the uncertainty on the source gives the major contribution. We 
could say that the transition between the small time regime and the large time 
regime occurs at time T when both contributions to L(t) are equally large, i.e. 
L(0) = A(pij(0) + p2j(0))T/m. Because we can at best reach the Heisenberg 
uncertainty, a lower bound for the time is given by T = L(0) 2 k c /c with k c the 
wavevector corresponding to the compton wavelength of the particles. In the 
next section we will plug some numbers into this relation to show that we are 
not making a trivial point in stressing this additional error contribution. Let it 
for the moment be sufficient to provide an example which illustrates the impor- 
tance of this uncertainty on the source. We can use a Gaussian distribution to 
represent the momentum correlation 

t-,/ \ (P1+P2) 2 , . 

F(pi,p 2 )~e - . (17) 
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The smaller the value of a, the better the momentum correlation between the 
two fragments. In the limit a — > this distribution approaches the Dirac delta 
distribution. Remark that this distribution is not peaked around a certain en- 
ergy Eq as should be required for a decaying system at rest. In the Appendix it 
is explained why we can leave this restriction on the energy aside without chang- 
ing the main result. It will follow that a reasonable energy width, peaked around 
E , will imply only a minor broadening of the wavefunction. The wavefunction 
at t = is 

1 Mr 1 ,r 2 ,0)~ ( 5(r 1 -r 2 ) e - rW4n2 . (18) 

Thus clearly ip represents a decaying system because initially ri = r 2 . But 
for small values of a (when F is peaked around Pi + p 2 = 0) the probability 
of finding the particles at t = at some configuration is totally smeared out, 
although the probability has a maximum at the origin of the coordinate system. 
Note that although the relation dS/dpi = is satisfied, this condition doesn't 
restrict the place of decay near the origin of the coordinate system, as was 
assumed by Peres. So in the short time regime, it may be hard to speak of 
possible opposite movements of particles relative to the origin because we cannot 
exactly say (at least without measurement) where the decay of the system took 
place. As time increases the deviation from the detection probability peak at 
mil - ! + m 2 r 2 = will even increase with time as was shown above. However, 
because the distance from the particles to the source increases, the uncertainty 
of the source will become less important for the angular deviation. In the long 
distance regime the angular deviation will then be dominated by the momentum 
uncertainty. 

In section 4 we will show that we can retain the classical picture of a decaying 
system in the previous example, in both the short time regime and the long 
time regime, when it is described by Bohmian mechanics. Namely in the case of 
the considered wavefunction, the Bohmian particles will depart near eachother 
and will move along opposite directions. However, the place of departure will 
vary from pair to pair over an extended region. The more this initial region 
is confined, the less perfect the momentum entanglement will be, and the less 
perfect the opposite movements of the Bohmian particles will be. 

3 Opposite momenta and opposite directions in 
the case of the experiments of Popper and 
Pittman et al. 

The importance of the error contribution arising from the uncertainty of the 
source becomes clear if we consider the example of Popper's experimental pro- 
posal. It can easily be seen that Popper's experiment cannot occur in the large 
time regime. If the particles would travel long distances (and hence obtain 
good angular correlation), the virtual slit (which is of the order of the transver- 
sal deviation) would be too large to have virtual diffraction of the particles on 
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Bob's side. Hence Popper's original experimental proposal should be consid- 
ered in the short time regime. However, in this regime the uncertainty on the 
source becomes the most important contribution to the angular deviation and 
as follows from the analysis of Collett and Loudon, this uncertainty makes a 
detectable virtual diffraction impossible. This also implies that in discussing 
Popper's experiment one should be carefull with statements such as ". . . the al- 
lowed deviation from perfect alignment is of the order of A\pi + P2\l\pi — P2I, 
which is much too small to be of any consequence in the present discussion." 
and "... nearly perfect alignment can be taken for granted, ..." [4]. 

In the case of optical imaging, performed by Pittman et al., Peres' analy- 
sis can be applied. Let us first briefly review this experiment. The experiment 
uses momentum correlated photons resulting from spontaneous parametric down 
conversion (SPDC). In this process a pump photon incident on a nonlinear beta 
barium borate (BBO) crystal leads to the creation of a signal and idler photon. 
The place of creation within the BBO crystal is unknown. Due to momentum 
conservation the sum of the momenta of these photons has to equal the mo- 
mentum of the pump photon. This results in the momentum entanglement of 
the two photons, because the momenta of the idler and signal photon can be 
combined in an infinite number of ways to equal the momentum of the pump 
photon. In the same way the energies of the created photons add up to the 
energy of the pump photon. In the experiment the signal and idler photons are 
sent in two different directions where coincidence records may be performed by 
two photon counting detectors. A convex lens, with focal length /, is placed 
in the signal beam in order to turn the momentum correlation of the created 
photons into spatial correlation (see Fig. 1). In front of the detector for the 
signal beam an aperture is placed at a distance S from the lens. By placing the 
detector for the idler beam at a distance S' from the lens, prescribed by the 
Gaussian thin lens equation, i.e. 

l + i = 7 (19) 

and scanning in the transverse plane of the idler beam, an image of this aper- 
ture is observed in the coincidence counts. This image obeys the classical lens 
equations in the following sense. If a classical pointlikc light source would be 
placed in the plane of the aperture, where the signal photon was detected, it 
would have an image where the idler photon was detected. This is the spatial 
correlation of the photons. 

The experiment can be seen to occur in the long distance regime because of 
the lens. In some sense the lens can be seen as projecting the angular correlation 
at infinity, to finite distances (at distances 2/ from the lens). The better the 
momentum correlation, the better the angular correlation is at infinity or the 
better the optical imaging is. 

It is interesting to use the data of this experiment to give an example where 
the border can be situated between the long distance and the short distance 
regime for a stongly momentum correlated system. A lower bound for the time 
is T — L(0) 2 k c /c, hence if there would be no lens in the Pittman experiment, 
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S=2f S'=2f 

Figure 1: Unfolded schematic of the experiment. Particle 1 (the signal photon) 
is moving to the right and particle 2 (the idler photon) is moving to the left. 
We have put S = S' = 2f. Similar schematics are shown in [8]. There the 
trajectories represent conceptual photon trajectories, here the trajectories are 
the Bohmian trajectories of massive particles. The place of creation of the 
photons within the BBO crystal should be taken sufficiently extended in order 
to preserve the opposite momenta correlation. 

the distance travelled by the photons at that time would be R — L(0) 2 k c ~ 70m 
(L(0) = 2mm is the width of the pump beam and the wavelength of the pump 
photon is 351. lnm). This means that if we would create momentum entangled 
photons via SPDC, then the angular deviation would be dominated by the error 
arising from the uncertainty of the source within a distance of R. 

4 Bohmian picture of a decaying system 

In this section we give the Bohmian description of a system consisting of two 
particles with opposite momenta, resulting from a decaying system at rest. 2 We 

2 Our simplified Bohmian approach to a decaying system is to be distinguished from the 
one studied by Y. Nogami, F.M. Toyama and W. van Dijk [14], where a decaying system is 
represented by a particle that leaks out from a region surrounded by a repulsive potential 
barrier. 
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use the wavefunction in Eq. (3) with momentum distribution 

F(pi,p 2 ) = N6( Pl +p 2 )e- a Pi/ h (20) 

where TV is a normalisation factor. The parameter a sets the scale of the initial 
separation of the two particles, as will be seen soon. A small a will correspond 
to the considered physical situation, i.e. a decaying system. The parameter 
is introduced in order to avoid singularities arising from the delta distribution 
when calculating the Bohmian trajectories later on. Again we have left aside 
the restriction on the total energy, as required in the case of a decaying system 
(see Appendix) . Remark that the exponential factor in Eq. (20) doesn't restrict 
the value of the total energy for a small a. The wavefunction corresponding to 
the distribution F is 

with \x the reduced mass of the fragments: ji — ^ + At t = the 

probability distribution is 

|^( ri ,r 2 ,0)| 2 = A 2 ^^V( ri - r2 ) 2 / 2 ^. (22) 

It follows that a small value for a corresponds to the considered physical situ- 
ation of a decaying system. However, the place of decay is unknown. This is a 
consequence of the Heisenberg uncertainty, as explained in the section 2. 

The Bohmian trajectories Rj(t) of the particles are found by solving the 
differential equations 



dR 3 _ 1 Reip* (ri , r 2 , t)p J 'ip(r 1 , r 2 , t) 



(23) 



dt rrij |V>(ri,r 2 ,i)l 2 _ ^ 

Because (pi + p 2 )V' = the trajectories of the particles satisfy 

^(miRi+m 2 R 2 ) = 0. (24) 

This shows that the particles have opposite speeds and thus move in opposite 
directions. Integration of the differential equations (23) leads to 

Ri(t) = Ci + C^W + a 2 

R 2 (i) = Ci -C 2V /£ 2 /V + a 2 (25) 

where Ci and C 2 are arbitrary constant vectors. It follows that the particles 
also move along straight lines. If we consider an ensemble of identically prepared 
systems, the probability density of the Bohmian particles equals the quantum 
mechanical distribution |?/>| 2 [10,11]. This is also called the quantum equilib- 
rium hypothesis [15] and ensures the empirical equivalence between Bohmian 
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mechanics and standard quantum mechanics. This equality determines the dis- 
tribution of the vectors Cj in the ensemble. Because the probability distribution 
at t — is sharply peaked at ri = r 2 (for small values of a) , the particles depart 
near each other. As follows from Eq. (25) their further propagation proceeds 
along straight lines, in the direction of their connecting line. Thus opposite 
momenta lead to opposite directions of movement for Bohmian particles. But 
their place of departure is located within an extended area, in order to preserve 
momentum correlation. 

By using Bohmian mechanics, we are thus able to retain part of the clas- 
sical picture of a decaying system at rest. Remark the similarity in language 
with the one used by Popper to describe his experiment. The difference is that 
Popper assumed the particles to depart from a confined region (which is how- 
ever incompatible with opposite momenta in quantum mechanics). Note that 
the opposite movements of Bohmian particles can not be verified experimen- 
tally. This is because a measurement of the place of departure of the Bohmian 
particles would cause a change of the wavefunction (this is the collapse of the 
wavefunction in standard quantum mechanics), which in turn implies different 
particle trajectories. 

5 The experiment of Pittman et al. revisited by 
Bohmian mechanics 

Although the experiment of Pittman et al. can be correctly explained with 
quantum optics, we will provide a Bohmian account of the experiment, when it 
is "translated" into its massive particle equivalent. One of the reasons to use 
Bohmian mechanics is that it justifies the conceptual photon trajectories drawn 
by Pittman et al. [8]. I.e. the photon trajectories coincide with the trajecto- 
ries of Bohmian particles in the massive particle equivalent of the experiment. 
This (Bohmian) quantum mechanical approach is contrary to the explanation 
in terms of "usual" geometrical optics used by Pittman et al. In quantum optics 
these paths are usually regarded as a visualisation of the different contributions 
to the detection probabilities. 

Because there is at present no satisfactory Bohmian particle interpretation 
for photons, see discussion elsewhere [16-19], we will follow Peres' departure 
and consider the nonrelativistic massive particle wavefunction in Eq. (3) in- 
stead of quantum optics, to study the experiment by Pittman et al. The SPDC 
source then corresponds to a decaying system at rest, resulting in two energy 
and momentum correlated fragments. We will assume the total momentum of 
the fragments to be zero, instead of some fixed value corresponding to the initial 
momentum of the total system (which would represent the momentum of the 
pump photon) . This assumption corresponds to the "unfolded" schematic intro- 
duced by Pittman et al. [8] (see Fig. 1). In this way we can use the momentum 
distribution F defined in previous section 

F(pi,p 2 ) =NS( Pl +p 2 )e~ ap i /h . (26) 
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In previous section we described the free evolution after the decay of the system. 
The unknown place of decay in the massive particle case corresponds to the 
unknown place of creation within the BBO crystal in the photon case (which 
is also related to the width of the pump beam). To complete the Bohmian 
description of the massive particle equivalent of optical imaging, we just have 
to describe the system's interaction with the lens. In classical optics we can 
use ray optics to describe the action of the lens on an impinging light beam 
[20] . The rays are such that the Gaussian thin lens equations are satisfied. Two 
generic examples, which we will need later on, are the following. The effect of 
the lens on a plane wave is to turn it into a converging wave, with focus in 
the focal plane, such that the corresponding rays obey the lens equations. The 
characteristics of the converging wave are then determined by the momentum of 
the incoming plane wave and the focal length. A second example is a spherical 
wave, representing a point source. If we assume that the light source is located 
in a plane at a distance S from the lens, then the spherical wave will turn into 
a converging wave with focus in the plane at a distance S' from the lens so that 
1/S + 1/S" = 1// and the source, the image and the centre of the lens will be 
aligned. In massive particle quantum physics the equivalent of optical lenses 
are electrostatic or magnetic lenses. These electromagnetic lenses are generally 
used to collimatc or focus beams of charged particles. This field of research is 
usually called optics of charged-particle beams or the theory of charged-particle 
beams through electromagnetic systems. Most of the literature deals with the 
classical description of the particles and only recently the quantum mechanical 
approach has been studied, see for example Hawkes and Kasper [21], and Khan 
and Jagannathan [22] and references therein. Here wc will not bother about 
the detailed analysis of particles passing through such electromagnetic lenses, 
and use directly, in the spirit of de Broglie, the analogy with classical optics. 
For example we can describe the action of an electromagnetic lens as turning a 
quantum mechanical plane wave into a Gaussian wave (we can take this as the 
analogue of the converging wave in classical optics, because a Gaussian wave is 
contracting before expanding), determined by the momentum of the incoming 
wave and the focal length. This analogy is very appealing because the rays 
in classical optics can be "identified" with the Bohmian trajectories. This is 
because in the one-particle case, the curves determined by the normals of the 
wavefronts of the quantum mechanical wavefunction are just the possible orbits 
of a particle described by Bohmian mechanics [23]. If wc apply this to our 
decaying system, then every plane wave of the particle impinging on the lens, 
say particle two, in the integral in Eq. (3) is turned into a particular Gaussian 
wave. The resulting wave is then 

V'(ri,r 2 ,i) = yF(p 1 ,p 2 )e J (P 1 - ri -rf t / 2mi )/ ;i G(r 2 ,p 2 ,/)dp 1 dp 2 (27) 

where G represents the Gaussian wave. This wave is guiding the particles after 
particle two passed the lens. To avoid unnecessary mathematical complications 
when calculating the Bohmian trajectories implied by the wave Eq. (27), we 
assume that the place of decay of the system is somewhere in the middle between 
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the lens and the detector on the right (where the idler photon arrives in the 
experiment of Pittman et. at). This corresponds to a BBO crystal placed in 
the middle instead of it placed near the lens, as in the experiment. When 
particle two arrives in the vicinity of the lens, particle one arrives in the vicinity 
of the detector placed on the right. We can describe the effect of the detector 
on the wavefunction as the first stage of a von Neumann measurement process 
(see for example Bohm and Hiley [24]). In this process the wavefunction gets 
entangled with a pointer which is represented by a wavefunction 4> with a small 
width. In this way the wavefunction in Eq. (3) turns into 

V"(n,r 2 ,y,i) = J F(p 1 ,p 2 )5(a-r 1 )</»(y,i,a)e i ( Pl -( a - r2 )- Bt )/ ft dp 1 dp 2 da 

(28) 

where we have written the factorization into position eigenfunctions 5(a — ri) 
of particle one explicitly. It is supposed that the interaction of the detecting 
apparatus with the decaying system lasted long enough in order to assure that 
the wavepackets 4>(y,a) are non overlapping for each a. During this process, 
the pointer particle has entered one of the packets <fi(y, a), determined by its 
initial position. While the packets are non overlapping, only the considered 
wavepacket is determining the subsequent trajectory of the pointer particle. 
The other wavepackets can then be dismissed for the further description of the 
system. In conventional quantum mechanics this process would be treated as 
a collapse of the wavefunction, which is the second stage of the von Neumann 
measurement process. As a result, the effective wave guiding particle two, is 
also reduced to the following superposition of plane waves 

V> 2 (r 2 ,i) = N J e i(p-( a -^)-P 2t / 2m 0- Q P 2 /R dp 

= N ( ^ y /2 e -(-r 2 )V4^ + ^)_ (29) 

\a + it/2m 2 J V ; 

The phase of this wave is 

So the wavefronts of the guiding wave of particle two are spheres with centre 
in a. Because the detectors are placed in planes at distances S and S' from 
the lens, with S and 5" obeying the Gaussian lens equation (19), this wave will 
result, after propagation through the lens, in a converging wave with focus in 
the plane at a distance S from the lens and where the focus is determined by the 
Gaussian thin lens equations. Hereby we used again the analogy with classical 
Gaussian optics. If for example S — S' and if the centre of the lens is taken 
as the origin of our coordinate system, then the focus will be at —a (see Fig. 
1). As a consequence of the quantum equilibrium hypothesis, particle two will 
be detected in the focus of the wave. Because we used a Gaussian to describe 
the converging wave, the Bohmian trajectories will not be straight lines, but 
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will be curved (for images see Holland [23]). The curvature will depend on the 
width in the focus of the Gaussian. In the limit of a zero width, however, the 
trajectories will approach straight lines, directed from the lens towards the focus 
of the wave. When the coincidence detections are considered, it will appear as 
if particle two departs from the place of detection of particle one. 

This completes the Bohmian analysis of the phenomenon of optical imaging. 
Before the fragments reach the lens, they move along straight lines from the 
place of decay. Note that this place of decay is not fixed, in order to guarantee 
the momentum correlation pi + p 2 = 0. When one of the particles reaches the 
lens, its direction of movement will change in accordance with the classical thin 
lcnse approximations. We assumed hereby the place of decay to be centered 
between the right detector and the lens. It can be expected, although it is not 
proven, that a random place of decay (for example near the lens) will lead to 
the same results in the Bohmian description of the experiment. 

6 Conclusions 

In conclusion we showed that Peres' analysis concerning the question to what 
extent opposite momenta lead to opposite directions, is only valid in the long 
distance regime. In the short time regime there is an additional source of an- 
gular deviation. On the other hand the statement "opposite momenta lead to 
opposite directions" is true in the Bohmian language. I.e. Bohmian particles 
travel in opposite directions when the wavefunction has eigenvalue zero for the 
total momentum operator (however from an unknown place of departure). We 
also showed that Bohmian trajectories can be used to gain insight into the 
phenomenon of ghost imaging as reported by Pittman et al. In fact Bohmian 
mechanics can be used to describe similar quantum optical experiments as well, 
such as Popper's experiment in the version of Kim and Shih [7] and the ex- 
periment reported by Strekalov et al. [25]. In particular, Bohmian mechanics 
justifies the use of the conceptual pictures of photon trajectories present in these 
papers. Note that the experiment of Kim and Shih is very illustrative for the 
need for perfect momentum correlation of the photons, or equivalently that there 
must be very little restriction on the place of creation of the photons to create 
a perfect image. This is because Kim and Shih failed in their original intention 
to perform Popper's gedankenexperiment, due to the restricted diameter of the 
pump beam used in the experiment. The imperfect momentum correlation then 
led to an imperfect optical image [6]. This is immediately obvious when we 
consider our Bohmian picture of optical imaging, because if the momentum cor- 
relation is imperfect, the Bohmian particles will not move in opposite directions 
before the system reaches the lens. 
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Appendix 



Conservation of energy requires that the energy of the total system equals the 
energy of the decaying system Eq. If this decaying system is initially at rest, 
Eq will be the rest energy of the system. Suppose now that we take a delta- 
distribution for this energy restriction i.e. 

F( Pl ,p 2 ) = f(p u p 2 )S(E - E ) (31) 

where / determines the momentum correlation (this is for example the distri- 
bution in Eq. (17) or Eq. (20)). If we take a distribution / which is real and 
symmetric, i.e. /(pi,P2) = /(— Pi,— P2) then the probability currents of the 
two particles are both zero for all times, i.e. 

mi 

This implies that the particles show no evolution. In Bohmian mechanics this 
represents particles that stand still, because the Bohmian speeds are defined as 
dRi/dt = ji/|?A| 2 - As a result the considered distribution in Eq. (31) doesn't 
correspond to a decaying system. We can resolve this problem by allowing a 
finite energy width centered around Eq. However, it will follow that, if the 
wavefunction displays strong momentum correlation in the sense that pi +P2 = 

0, the restriction to a small energy width only involves a minor broadening of 
the wavefunction, which implies that we can leave the restriction on the total 
energy aside for our qualitative analysis. 

For the momentum distribution / we will take the distribution in Eq. (20), 

1. e. /(pi,P2) ~ <5(Pi + P2)e~ QPl//n . The restriction on the total energy is ac- 
complished by integrating over values for (p 1; p 2 ) for wich E_ < E < E + , for a 
certain minimum energy value E_ and a certain maximimum energy value E + . 
We therefore define the following function 

j. / rp W * I 1 if + ^ < £± 

disc(.E±)(Pi> P2) = < 2m i 2m 2 - 1 

I otherwise 

The momentum distribution then becomes 

F(pi, P2 ) = /(pi, P2 )[disc(£ + ) -disc(£7_)]. (33) 

The resulting wavefunction of the system is then 

V(ri,r 2 ,£) - y /( Pl ,p 2 )[disc(ii; + ) -disc(S_)] e 4 (P 1 - ri +P 2 - r2 - £t )/ R dp 1 dp2 
~ y e ip - (ri - r2)/ft - (<t/2 " +a)p2/ft [disc'(-E+)-disc'(£;_)]dp (34) 

where 

1 if p 2 /2n < E± 



disc'(£±)(p) 



otherwise 
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If we write Eq. (34) as a Fourier transform then we can apply the convolution 
theorem 

^( ri ,r 2 ,t) ~ ^-r,)/n}^- {U/2 " +a)p2 ' h ) ® ^ ri - r2 )/ fi }( disc '(^ + ) - disc'(i? 
<~ h(x, t) <g> g(x) 

where we have defined 

MM) = ^ + (ri _ r2)/ , } ( e - (jt/2 ^ )p2/?l ), 
g(x) = ^ { + (ri -r 2 )M}( disc '^+)- disc '^-))' 

x = |ri - r 2 |. (36) 

We present now two methods to show that the restriction on the energy can 
be relinquished, without changing the qualitative analysis. The first method 
proceeds as follows. The first function in the convolution in Eq. (35) is just the 
wavefunction in Eq. (21), 

If we define a± — 27r^/ E±2fi/h then the second function in the convolution in 
Eq. (35) in two-dimensional physical space becomes 

g(x) ~ [a + Ji(a + x) — a_ Ji(a_x)] jx (38) 

where J\ is the first order spherical Bcsscl function. To evaluate g(x) we will 
substitute some reasonable values for E+ and E- in Eq. (38). In addition we 
will assume that the fragments have equal masses so that we can put 2[i = m, 
with m the mass of one fragment. For E+ we will take one percent of the rest 
mass of the total system in order to avoid the relativistic regime, E + = 0.02mc 2 . 
We will take an energy gap of 0.0011S+, so that £7_ = 0.999-E+. In this way 
a + « 5.58309/A c and a_ « 5.58023/A c , where A c is the Compton wavelength 
of the fragments. In Fig. 2 the function g(x) = \a+Ji(a+x) — a_ Ji(a_x)] jx is 
plotted for x in units of the Compton wavelength A c . 

The figure shows that g(x) is a rapidly oscillating function with a peak at 
x = 0. We now give the distribution of h(x,t) at t — a width of the order 
of twenty times the Compton wavelength, which can be done by adjusting a. 
Recall that the function h was in fact the wavefunction of the system if we 
didn't restrict the energy. So the width of h is in fact the measure of the 
initial nearness of the fragments, which is then of the order of twenty times the 
Compton wavelength. Then due to the rapid oscillation, the main contribution 
in the convolution will arise only from the peak in g(x) at x — 0. This peak will 
result in only a small broadening of h(x, 0) so that h(x, 0) ® g(x) w h(x, 0). 

Because the width of h(x, t) only increases with time, this approximation 
will become more precise with time. In conclusion we can put the unnormalised 
wavefunction equal to 
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x in units of the Compton wavelength 



Figure 2: The function g(x) = [5.58309 Ji(5. 58309a;) - 5.58023 Ji (5.58023a;)] /x 
is plotted for x in units of the Compton wavelength A c . 



The larger the energy gap or the larger the rest energy, the more rapid the 
oscillation of g(x) will be and the more peaked g(x) will be at x = 0, and as a 
result the more the approximation is valid. 

A second way to achieve this result is to assume that h(x, t) is very narrowly 
peaked at x = for t — 0, so that h(x, 0) « S(x). This is the case if a approaches 
zero. As a result h(x, 0) ® g{x) « g{x). Thus in this case the unnormalised 
initial probability distribution is g{x) 2 . This distribution is plotted in Fig. 3 
with again x in units of the Compton wavelength. This figure shows that most 
of the probability is concentrated within a few times the Compton wavelength. 

This implies that the particles depart from a very narrow region, only a few 
times the Compton wavelength in diameter, from each other. However, it is a 
priori unknown wether the Bohmian particles will travel along straight lines, 
in the direction of their connecting line, in this case . But this can be easily 
proven by invoking symmetry arguments. Because the wavefunction ip(ri, r2, t) 
in Eq. (35) is only dependent on the difference |ri — r 2 |, the Bohmian velocities 
(as defined in Eq. (23)) will be opposite and the Bohmian particles will travel 
along straight lines in opposite directions. 

In conclusion, we have shown that the energy restriction does not put a 



18 



0.0003 



0.00025 



0.0002 



0.00015 



0.0001 



5e-05 




x in units of the Compton wavelength 



_. „ ml . . . [5.58309Ji(5.58309a;)-5.58023Ji(5.58023a;)l 2 . , , 

figure 3: Ihe function g(x) = - ^ — is plotted 

for x in units of the Compton wavelength A c . 



restriction on the Bohmian picture of the system. The only effect of the energy 
restriction is a minor broadening of the probability density. 
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